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In the zero temperature Glauber dynamics of the ferromagnetic Ising or g-state Potts model, the 
size of domains is known to grow like t^''^. Recent simulations have shown that the fraction r{q,t) 
of spins which have never flipped up to time t decays like a power law r{q,t) ^ i"*'?' with a non- 
trivial dependence of the exponent 6{q) on q and on space dimension. By mapping the problem on 
an exactly soluble one-species coagulation model {A + A ^ A), we obtain the exact expression of 
6{q) in dimension one. 
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Phase ordering and domain growth in systems quenched from a disordered phase to an ordered phase has been a 
subject of much interest during the last fifteen years in fields ranging from metallurgy to cosmology Q. It is well 
established that the pattern of growing domains is self-similar in time and that the characteristic domain size increases 
with a simple power law f. For example, p ~ 1/2 holds for all systems with short range-interactions described by a 
scalar non-conserved order parameter. However, as noted for the auto-correlation function correlations at different 
times are characterized by more complicated exponents. Recently, the fraction of spins r(q, t) which have never flipped 
up to time t has been measured in simulations of coarsening at zero-temperature for the Ising and for the q-state 
^ Potts models. This fraction decreases with time like a power law [00], 

r{q,t)r.t-^i^) (1) 

o ■ 

, and numerical data indicate that the exponent 9{q) varies both with q and the dimension of space. Since dr{q,t) / dt 
' measures the probability that a given point is crossed for the first time at time t by a domain wall, 9{q) can be viewed 
'n . as a first passage exponent 

The aim of the present letter is to give the exact expression of 6{q) in one dimension. 
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This result fully agrees with previous numerical predictions based on MonteCarlo simulations ||3|,^ or on finite size 
scaling calculations It implies that for the Ising model 9{2) = 3/8 is exact. Note however that for other choices 
k>( ' of q, the exponent 9{q) is in general irrational (for example 0(3) = .53795082..). 

; I To obtain (H), we are going to follow four main steps: first by using finite size scaling we will relate the exponent 

\ 9{q) to the large L behavior of the fraction pl{<i) of spins which never flip between time and time oo for a finite 
one dimensional system of L sites with periodic boundary conditions (for the zero temperature Glauber dynamics of 
the g-state Potts model); secondly, we will show that the calculation of pl((7) can be reduced to solving the steady 
state of a reaction-diffusion model {A + A ^ A) on a one dimensional lattice of L sites with a source of particles at 
the origin (i.e. at site = L); our third step will be the solution of the steady state of that reaction-diffusion model 
leading to the exact expression of Pl{q) for arbitrary L and q; lastly, we will extract the exponent 9{q) from the large 
L behavior of this expression. 

In an infinite system, it is known |^,^ that starting with a random initial condition, the size of domains grows like 
t^^^. For a finite system of size L, one expects the dynamics to be very much the same as for the infinite system when 
t ^ (as the size of domains is small compared to the system size). On the other hand, for t ^ L^, there is a single 
domain left in the system and the dynamics stops. Therefore one expects that pL{q) ^ r{q, L^) which implies that 
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We obtain, in what follows, an exact expression of pL{q) valid for all sizes L and (||) will allow us to extract then the 
exponent 9{q). 

We now show that the calculation of pL (q) is equivalent to solving the following reaction-diffusion ||] model (A-l- A 
A) defined on a ring of L sites with a source of particles at the origin (that we choose to be at site L = 0): the origin 
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(site L) is always occupied and the other sites i (for 1 < i < L — 1) are either occupied by a particle A or empty; during 
every infinitesimal time interval dt, each particle hops with probability dt to its right neighbor and with probability 
dt to its left neighbor (and does not move with probability 1 — 2dt); if two particles occupy the same site, they 
instantaneously coagulate {A + A A); in addition whenever the particle at the origin jumps to one of its neighbors, 
a new particle is instantaneously produced at the origin. 

In the steady state, injection of particles at the origin is balanced by aggregation of particles in the bulk and one 
has a probability Pi(r7i) of finding m sites occupied on the ring of L sites. The connection with the spin problem 
described above is made by expressing the probability pl{<i) of never flipping (from t = until t = oo) in terms of 
PL{m) 

PL{q) = E Pl{^) . (4) 

m=l ^ 

This formula is obtained by remembering that updating a spin in one dimension with Glauber dynamics at zero 
temperature simply consists in choosing for its new value the value of one of its two neighbors at random. When the 
value of a spin Si at time t is traced back in time, a random walk is obtained which connects Si{t) through various 
ancestors to a particular spin in the initial configuration with value V . Now, let us consider all the updates of the spin 
at the origin. They all give rise to random walks going backward in time. These walks can merge and are created at 
the origin exactly as in the one species coagulation model. So, in the limit t — > oo, they lead to m different ancestors 
in the initial configuration with probability PL{m). Equation (^) then follows by noting that the spins in the initial 
configuration are uncorrelated and that the spin at the origin never flips if and only if all its updates have the same 
value. 

We now come to the full exact solution of the reaction-diffusion model steady state which is the crucial point of 
the present work. It is known that some properties of the coagulation model A + A ^ A can be calculated exactly 
[0 even in the presence of a fixed source ||n[]. The properties which have been calculated jl^ so far on this sort 
of problems are the probabilities that a (connected) region of consecutive sites are all empty. This is because one 
can obtain closed kinetic equations (see (^ below) for these quantities. The key to our exact solution is that more 
complicated quantities like the probability of having 2 (or 3 or • • • n) disconnected empty regions can be expressed in 
a simple way in terms of these probabilities that a single connected region is empty. 

Let us define for the coagulation model the probabilities Bij, 1 < i < j < L that the segment {i, i + 1, j — 1} 
contains no particle . These quantities change when a particle enters or leaves the segment through its extremities 
and this can be expressed in terms of the Bi j themselves (for example, the probability to find a particle at site j 
given that sites i, z + 1, • • • , j — 1 are empty is Bij — Bi^^i). Therefore, the Sjj's satisfy closed kinetic equations 
which read in the steady state 

Bi+ij + Bi^ij + + Bij^i — iBij = (5) 

Modified equations for j = i + 1 and boundary conditions coming from the permanent occupation of the origin can 
be taken into account by the conventions Bij — 1 and Bq j — Bij^j^x = 0. 
The explicit solution of (||) is ||]: 



k even k' odd 



Ssinfca sinfe'a (sinfeia sinfc'jo! — sinfeja sinfe'ia) 
(L + 1)2 (cosfc'a — cosfca) (2 — cos fea — cos fc'a) 



with a — Tr/{L + 1), 2 < k < L, 1 < k' < L (Note that (||) is not valid for j — i since it gives Bu ~ instead of 1). 
For example, when i = 4, this gives Bi,2 = B^a = 26/44; Bi^a = ^2,4 = 16/44; Bi,4 = 8/44; ^2,3 = 30/44. 

These known results can be generalized by introducing the probabilities ... ^ that there is no particle 

in any of the disconnected segments + 1, ■ ■ ■ ,i2 — 1}, ■ ' ' {i2n-i, ■ ■ ■ , *2n — 1} with ii < ?2 < • ■ ■ < *2n- Similarly 

to the i?ij's (to which they reduce for n — 1), the B^^^'s satisfy equations analogous to (^). The new boundary 
conditions can be taken care of by the convention that S'") reduces to the corresponding _b("^i) if two successive 
indices coincide. 

It turns out that the diffusion- like equations satisfied by the i3("^'s can be solved exactly and that their solutions 
can all be expressed in terms of the Btj. Namely, 



BS,k,i = B,,jBk,i + B^,iBj,k - B^,kBj,i (7) 



and more generally 
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where the sum runs over all the permutations of the indices {ii, 12, ■ • ■ , i2n} and e{a) is the signature of the permutation 
a (Note that (0) and (||) are Pfaffians fl^). One way to prove is simply to check that their r.h.s. satisfy the 

same diffusion-like equations as the ^("'^when Bij is solution of (g|) and that they indeed reduce to the expression 
for if two successive indices coincide. In the case L — -i, one finds that b[\^ 4 = 15/44. 

Once the B*^"-' 's are known, all the steady state properties of the reaction-diffusion model are in principle computable, 
in particular the P£,(m), and as a consequence (^) the There is a simple way of obtaining the Pl{(i) in terms 

of the i3(")'s: introduce the (normalized) weights w{ti, ■ ■ - tl^i) in the steady state of configuration {ti, • • - tl^i} 
(where = 1 if site i is occupied and = if it is emtpy); we can write the sum over possible particle numbers in 
the steady state in eq.(^ as an explicit sum over all possible configurations 

yp,(,,) 1 = V "^"^v";""-^^ (9) 

m=l tx,---tl-i 

Using the identity q^-'ri+---+^-'rL-i — nj=i^[l + ^ 1)(1 ~ T?)]' expanding the product and regrouping the different 
terms, we obtain 
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l<i<j<L l<i<j<k<l<L 
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This gives the pL{q) in terms of the known Bij and it can be used in this form. For example, for L = 4, it gives 
Pa{<i) = (Sg^ + 23(7^ + 12(7 + l)/44g^ in agreement with what had been obtained by a direct calculation for small sizes 

An equivalent but more compact and manageable form can be obtained by defining two antisymmetric matrices 
Aij = —Aji = —{q — l)^^*and = — Cji = Bj.j for 1 < i < j < L. Eq.(p^ can then be rewritten simply as p^ , 

P^il) = ^Vdet(l + AC) (11) 

This expression ( pi] ) is our main exact result. 

We have only been able to compute exactly the determinant of eq.([Tl|) for q=0 (and l!) For other values of 
we could only evaluate its asymptotic behavior. We summarize the main elements of our estimation which will 
be detailed elsewhere jl^ . For large L, and i,j far from the boundaries, Bi j becomes a continuous function (3{x,y) 
oi X = i/L and y = j/L which satisfies the Laplace equation (the continuous version of eq.(|^)) with the boundary 
conditions /3(x, x) = 1, /3(0, x) = P{x, 1) = 0. This is a slowly varying function except in the two symmetric arbitrary 
small corners < x < ?; ^ 1 and 1 — j/<l — x^l where its asympotic behavior is 

l3{x,y) ^ (3{1 ~ y,l ~ x) = -ta.n^^{x/y), for < x < y < 1 (12) 

TT 

These two small neighbourhoods are responsible for the singular behavior of the determinant for large L. In the region 
a; < y <C 1, (3{x,y) is of the form f{x/y) and the singular contribution can be obtained by computing the traces of 
the power of a Toeplitz matrix in the variables \og{x) ,\og{y) . The same is true in the region 1 — j/< 1 — x 
where /3(x, y) ~ /[(I — 2/)/(l — x)] (and with the variables log(l — x), log(l — y)). In this way we obtain for 1 < q <2, 



0{q) = I £ log (1 - 2^ dx {x^'^ + x-^") fix) ] (13) 







Using (p^), this leads to (|^) for 9(q). When q > 2, there is an additional term to ( p^ coming from an eigenvalue of 
the matrix 1 -I- AC which becomes very small as L — > 00. As a consequence, ( p^ has to be replaced by its analytic 
continuation from the range 1 < g < 2 and this leads again to @). 

In the whole range of q, the result (^) is in excellent agreement with previous numerical estimates For 
example, for q = 5, the exact prediction (^ gives 6{5) ~ .6928365.. which agrees with the numerical estimate 
.6928 ± .0003 of 0. 
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Our result (|^) shows that systems as simple as the one dimensional Ising or Potts model at zero temperature 
exhibit rather complicated power law decays very similar to other irreversible processes such as reaction-diffusion 
or sequential-parking Jl6t problems. 

Result (^) is also reminiscent of the exactly known exponents jl^ ] of the equilibrium Potts model in dimension 2. 
It would be interesting to see whether the methods originally used for these equilibrium problems or the more recent 
conformal theory techniques could be extended here to rederive (|^). 

We have already noted that expressions (0) , (||) , (|l^) and (|l^) are Pfaffians which are always present in free 
fermion problems. In fact, another way of solving the coagulation model (with its boundary) consists in mapping it 
onto a free fermions problem. The details of this solution which requires several successive transformations of the 
original problem will be given in As the problem of the steady state of the reaction-diffusion is completely solved, 
one could calculate steady state properties such as density profiles, correlation functions. Also, as the whole problem 
can be reduced to a free fermion problem [Q (through unfortunately rather complicated transformations) , one should 
be able to calculate all kinds of unequal time correlations. 

Lastly, it is worth noting that the different scaling behaviors of pL{q) when q varies, arise in (^) from the asymptotic 
form of PL{m) for L ^ 1, PL{m) ^ exp[log-L g{m/ logL)] where the function 5 is a Legendre transform of 2d{q) with 
respect to log(<7) in a way quite reminiscent of the usual description of multifractal measures |l8|| . 
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